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Abstract 

The Casimir energy for a compact dielectric sphere is considered 
in a novel way, using the quantum statistical method introduced by 
H0ye-Stell and others. Dilute media are assumed. It turns out that 
this method is a very powerful one: we are actually able to derive 
an expression for the Casimir energy that contains also the negative 
part resulting from the attractive van der Waals forces between the 
molecules. It is precisely this part of the Casimir energy that has 
turned out to be so difficult to extract from the formalism when using 
the conventional field theoretical methods for a continuous medium. 
Assuming a frequency cutoff, our results are in agreement with those 
recently obtained by G. Barton [J. Phys. A: Math. Gen. 32, 525 
(1999)]. 

KEY WORDS: Casimir energy; van der Waals forces; Quantum statistical 
mechanics; Polarizable fluids; Radiating dipole interaction. 

^^Written for a festschrift issue of Journal of Statistical Physics , dedicated to George 
Stell. 
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1 Introduction 



It is a pleasure to contribute this article to a festschrift volume for Professor 
Stell. The article is one in a series of articles published by the present au- 
thors on the Casimir effect and related topics, using methods of statistical 
mechanics for quantized systems at thermal equilibrium. Besides contribu- 
tions from others, these methods were developed by one of the authors ( JSH) 
in cooperation with Professor Stell in their extensive studies of polar and po- 
larizable fluids through several years. This long-lasting cooperation, which 
was initiated in 1972, is still active today. 

The Casimir energy problem for a compact spherical ball is a many- 
facetted problem; the formal solution of it is to an unusual degree dependent 
on the mathematical method of approach chosen. The Casimir effect as such 
is now a well-known effect in physics [p]]. It is ordinarily examined with the 
use of field theory in dielectric media, allowing the medium to possess a re- 
fractive index n (even dispersive effects can in principle be dealt with in this 
way, if n is assumed to depend on the frequency). The standard configura- 
tion does not involve curved boundaries at all, but consists instead of two 
plane plates separated by a small gap. In this geometrical configuration the 
phenomenological electromagnetic theory, as constructed mainly by Lifshitz 
[0] is fully adequate, and leads to a prediction for the Casimir force between 
the plates that has recently been verified experimentally to an impressive 
accuracy of about one per cent 0, ||. 

If we now leave the parallel-plate configuration and consider instead a 
single dielectric ball, the situation becomes much less clear-cut. The history 
along this direction of research may be taken to start with the calculation 
of Boyer on a singular perfectly conducting shell ||: he found the Casimir 
energy E to be positive, corresponding to an outward directed surface force. 
Later, the dielectric ball was considered by Milton ||, Milton and Ng 

[Pl,fll0|, Brevik et al. fni,[]P2"|,|r5|l,[]l}| and several others. Some consensus 



seems by now to have been reached as regards the Casimir energy E as 
found by field theoretical methods: this energy is positive, corresponding to 
a repulsive force, and is given by 

E= 23_hc ( 1)2 (1) 
3847ra V ; ' V ; 

for a dilute sphere whose radius is a. 
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Faced with this field theoretical result one becomes however surprised, for 
the following physical reason: the Casimir energy should be the cooperative 
result of the van der Waals forces between the molecules in the ball. The 
van der Waals forces are necessarily attractive. How can these forces sum up 
to give a repulsive total surface force? The natural answer to this question 
is that the field theoretical calculation, based as it is on a continuum model 
for the dielectric, is unable to copy with the attractive part. In other words, 
the attractive terms are necessarily lost in the regularization process. An 
important progress was recently made by Barton [pL5| ; he made use of quan- 
tum mechanical perturbation theory to second order, imposed an exponential 
cutoff in wave numbers, and arrived at a definite expression for the Casimir 
energy containing also the cutoff dependent, attractive (and actually also re- 
pulsive) terms. Moreover, a cutoff independent, repulsive term was contained 
in the energy expression, which was in precise agreement with Eq.(l) above. 
There are actually some indications of the same kind already in the paper 
of Milton and Ng ||10|| : they derived the cutoff independent Casimir energy 
starting from the van der Waals forces, omitting the divergent terms. 

And this brings us to the central theme of the present paper, namely to 
rederive the expression for the Casimir energy using the perhaps somewhat 
more unconventional quantum statistical methods that were developed by 
H0ye and Stell, and others. Central references for the present work are ||16| 



and fL7| . Others that also include evaluation of frequency spectra are |1S 



As we will see, this method is very powerful, and we will be able to establish 
contact with the results of Barton. The line of development of the application 
of this method to the Casimir problem is the following: Some years ago 
Brevik and H0ye |19] showed that the Casimir energy between two point 
particles is the same as the free energy due to two quantized fluctuating dipole 
moments interacting via the dipolar radiation interaction (zero frequency 
limit corresponds to the static dipole - dipole interaction). Later, H0ye and 
Brevik ]2D| extended this method to evaluate the Casimir force between a 
pair of parallel dielectric plates separated by a small gap. Performing this 
more complex calculation with the use of statistical mechanics for systems 
in thermal equilibrium, we were able to rederive the known results. 

Below we will evaluate the free energy in a dielectric dilute medium, again 
using the same methods of statistical mechanics. Based on this we will make 
contact with the results of Barton, as mentioned, as well as with the results 
obtained in field theory. On the basis of our method the physical origin of the 
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divergences is easily understood. The problem, as anticipated above, has its 
origin in a continuum description of dielectric media, while a realistic system 
has to have a microscopic structure involving a minimum separation between 
molecules due to repulsive cores. 



2 Basic formalism 

We begin by recapitulating some of the basic formulas from our earlier work 



HJ . For a pair of polarizable particles the free energy F due to their mutual 
interaction is 

-/^ = |E«x(2^(r)+^(r)); (2) 

Z K 



cf. Eq.(5.14) in |T9|]. Here «jf is the frequency dependent polarizability, and 

K = 2im/p (3) 

is the Matsubara frequency related to the frequency u via 

K = -ihuj. (4) 

Further, n is an integer, (3 = l/ksT is the inverse temperature, and i/jdk and 
ipAK are the two radial parts of the radiating dipole - dipole interaction as 
given by Eqs.(5.9) and (5.10) in fL9| . Performing the sum in (2) for f3 — > oo 
(i.e. T -> 0), we obtain Eq.(5.16) in Q 



23hca 2 

(5) 



47rr 7 

which is the known result for the Casimir effect. 

For a low density medium the total free energy A can now be obtained 
by summing or integrating (2) over pairs of particles in a volume V such that 

A = ip 2 |dr lC /r 2 F, (6) 

where p is the number density and r = ri — T2- At T = 0, (5) is to be inserted. 
Clearly, the integral will diverge, due to the behaviour for small values of r. 
However by integrating over a small sphere a finite term, which is positive, 
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can be separated out; cf. This finite term turns out to coincide with the 



field theoretical result. We will show below that the divergences found using 
other kinds of approach are connected with this small r behaviour. Equation 
(2), together with (6), will be shown to lead to results in agreement with those 



obtained from quantum mechanical perturbation theory to second order [|15 



The exponential cutoff used by Barton in Fourier space can be introduced 
also in our approach; it corresponds to a small "soft" r cut out from the 
otherwise continuous medium, and will be a rough approximation to real 
systems. As mentioned above, real systems are not continuous but consist of 
molecules that have a minimum separation due to hard cores. 



3 Calculation of the free energy 

Let us now calculate the free energy A, as given by (6). In order to establish 
connection with perturbation theory, we first represent (2) in terms of Fourier 
quantities for which, as we will see, a wave vector cutoff can be introduced 
easily. The radiating dipole - dipole interaction used in (2) can be written as 

V(12) = ^ DK {r)D K {\2) + VA*-(r)Atf(12), (7) 

with 

D K (12) = 3(ra 1K )(ra 2K ) - ai K a 2K , 
A^(12) = a 1K a 2K . 

Here the hats denote unit vectors, and dux is the Fourier transform of the fluc- 
tuating dipole moment of particle number % in imaginary time; cf. Eq.(5.2) 
in fl9fl . Equation (7) can be Fourier transformed to give 

$(12) = 4>AK{k)D K {l2) + 4> AK {k)A K {12), (8) 

with 

D K (12) = 3(ka 1K )(ka 2K ) - a 1K a 2K , 
V>(12) = 7335/ $(12)e lkr rfk. 
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With this Eq.(2) can be rewritten as 



-/3F=^E/M,e<««> dk *' (9) 

where 

M K = $ AK (k)ip AK (k')(3(kk>) - l)+i> AK (k)^ AK (k'). (10) 

Like expression (2), this is obtained after orientational averaging of the 
products of terms containing Dk(12) and A#(12) with respect to am- 

To get further the explicit Fourier transformed interaction terms are 
needed. These follow from the solution of Maxwell's equations, and like 
the corresponding terms in (2) (Eq.(5.10) in |19|]) they were used by H0ye 



and Stell when dealing with the refractive index of fluids [T7|. Thus from 
Eq.(7) in [P7] we have [] 

~ „x 4vr k 2 ~ 4?r / 2k 2 A 
fe(fc) = -^tt^ ? , Mk) = — j-, ? - (2 + 9) , (11) 



3 A; 2 - 

with = —ihcuj. (For simplicity is replaced by cu; where c is the light 
velocity.) Here the parameter 9 introduced by H0ye and Stell [^J is used. A 



purpose to introduce it was by 7-parametrization of the dipole-dipole inter- 
action to obtain a continuous family of mean field theories (7 — - > 0) of polar 
fluids. Here 7 is the inverse range of i/j A k(t), and for i(jdk{t) it is the inverse 
range inside which the dipolar 1/r 3 behaviour is cut or rounded off. As seen 
from (11) the 9 is thus the integrated amplitude of the '(/'Ait-term (u = 0). 
This parameter was also used in [IT]], part IV, and its Eq.(57) for the di- 



rect correlation function corresponds to Eq.(ll) here. With its Eqs.(56) and 
(59) the dielectric constant e can then in general be expressed as (for small 
w/7 -> 0) 

e-l Ait 
(l-9),+ (2 + 9) = Y pa > (12) 

where p is the number density of particles. (Here a possible density depen- 
dence of a — > a e ff which is proportional to the fluctuating dipole moment 
squared, is disregarded.) 

2 Note here that the Fourier transform in imaginary time on the interval from to (3 is 
the same function as the real time transform; cf. the derivation in Appendix B of jl{§. 
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The separate term (2 + 0) at the end of (11) will necessarily yield infinity 
when inserted in (9) and summed. In r-space it gives a 5-function at r = 0. 
As pairs of particles in reality are separated, we can simply remove this term 
here. This amounts to putting 2 + O = 0, by which 
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An pa. 



(13) 



Note here that this choice is consistent with the continuum approach (7 — > 
00), where only a transverse radiating field is implicitly considered. That 
is, with O = —2 the longitudinal part vanishes, as follows from Eq.(28) 
in DJ. Then one has (with the direct correlation function c — > if>) ipi — 
I^ak+^dk) = 0, while the transverse part becomes ip2 = \ijpAK~ 4 s dk) = 
(An/3)k 2 /(k 2 -lu 2 ). 



Inserting (11) into (10) we obtain 
(An) 2 k 2 k' 2 



M, 



K 



3 k 2 — uj 2 k' 2 — uj 2 
(An) 2 k 2 k' 2 r 1 



[(jfejfe 



'\2 



3 k 2 - k' 2 



+ 1] 
1 



k' 2 



k 2 



(kk'Y + i 



(14) 



Now the summation in (9) can be considered, and by restricting ourselves 
to frequency independent polarizability ax = ot we can easily perform the 
sum since (14) is of standard form for simple harmonic oscillators |16[ . We 
have (K = 2nn//3) 



x - (M)) 2 
^ (^ ) 2 + K 2 



1 cosh(|/37ia;o) 
r hU ° sinh(i/?^ ) 



(15) 



With ax = a, use of (15), and (14) inserted for M K , we now easily find at 
T = (K = —ihcu, hujQ — > hck and Tick') 



K 



(An) 



3 k 2 - k' 2 2 



(3tic(k'k 2 -kk' 2 )[(kk') 2 + l] 



(16) 
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When inserting this into (9), and further inserting into (6), we see that the 
result will diverge due to the small r (or large k) behaviour. This divergence 
can be avoided by introducing a large wave number cutoff, as Barton did [|T5|]. 
Thus we incorporate a factor exp(— Xk) in the interaction terms (11), which 
implies a factor exp(— X(k + k')) in (16). Regarding the electromagnetic field 
as a set of harmonic oscillators that mediate the interaction between the 
particles, the effect of this cutoff is to remove the high frequency oscillators. 
With (16) and (9) inserted into (6) we then obtain 

A = ~^y 2 J dr ld r 2 J dkdk'kk'e^'> ^^[(* k'f + 1], (17) 

where 7 is given by (13). 

With A = 7 2 A2 this is precisely the result obtained by Barton [15| ap- 



plying quantum mechanical perturbation theory to the dielectric continuum. 
That is, we have recovered the second order contribution which in the case 



of a scalar field is given by Eq.(4.2) in [15 where the equality 
exp(-X(k + k')) 



k + k 1 



/ dtex P (-Z(k + k')) (18) 



has been used. When dealing with the electromagnetic field, as we do here, 
the result in |15 is modified with the factor [(k k ) 2 + 1] as given by Eq.(B.2) 
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15] . With this, it is seen that our results are identical with those of 
Barton. Thus, further considerations based upon this basic agreement will 
necessarily be the same, and will not be repeated here. 



4 Further remarks 

Let us make a few remarks on the resulting free energy (or internal energy at 
T = 0) for a spherical body of radius a where a positive cutoff independent 
term going like 1/a shows up also within the field theoretical approach. This 
is commonly interpreted as equivalent to a repulsive Casimir surface force. 
From our approach it is now obvious that this term reflects the large r- 
behaviour of the free energy (5) of particle pairs. For a finite system such 
as a sphere the resulting free energy will in general be larger than its bulk 
value since there is no material outside with which it can interact. This 
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missing interaction first of all manifests itself in a positive term that reflects 
the surface tension and is proportional to the surface area. This term is 



present in the final result given in the Abstract in [TE\. In addition, there is 
the 1/a term which can be associated with the l/r 7 -tail in connection with 
the missing material away from the surface. 



As another point, let us note that in |15| there is an additional leading 



term due to first order perturbation theory. This term is not present in our 
derivation above. However it can be identified in a straightforward manner 
from the earlier work of H0ye and Stell fll7| . It represents a self-energy of 



the electromagnetic field attached to the polarizable particles. As such it is 
just part of the properties of isolated single particles, and should accordingly 
not be included in the free energy density above. In fact, it is part of the 
radiation reaction of accelerated charges and is incorporated in the resulting 
physical momentum. Unfortunately this momentum correction for a classical 
particle is infinite , so the "bare" mass of particles is negative resulting in 



the well known "runaway" problem in classical electromagnet ism ||22|| . This 



also reflects itself in the refractive index problem considered in jr?) . Further, 
H0ye and Lomba made numerical calculations to obtain a minor non-causal 
tail in the dielectric response of the fluid ( i. e., a minor response will appear 
before an electric field is applied to avoid exponential growth or runaway for 
increasing time p3|). 

With cutoff in wave vectors the above mass correction can be made finite. 
The term of interest is then the component ipAK{k) of Eq.(8). Including the 
shielding factor in (11) we then have (with 2 + © = as above) 

4tt 2k 2 
~3k 2 -u 2 



$AK{k) = —- ,e' Xk . (19) 



As shown by Eqs.(12)-(15) of fTTH , this gives rise to a self-interaction (for 
given K, s — > clk) 

A0(a x ) = -I4VWO), (20) 

where the minus sign is due to the definition used. The total internal energy 
contribution that follows from this is 

A 1 = 5>F<A0K)> = -IpVJ2(^ak(0), (21) 

K £ K 

where pV is the number of particles in a volume V . 
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Now (a 2 K ) = 3ak//3; cf. fl9[| . So with a;# = a we obtain by first using 
(15) for (3 -> oo 



or 

3 (2tt) 

which together with (13) and (21) yield 



E</W(0) = — 7^ /A;e A *dk=-^ r , 



1 Tr 3aAf3hc 



"7 



3 he, 
2^2 



V. 



(22) 



(23) 



This is precisely the first order contribution obtained by Barton ||15|| . As 
mentioned above, this energy is a quantity that is part of the free particles 
themselves, and can thus not be separated out. So the free energy of the 
interacting system of polarizable particles will not include this term. 

After the above was written, we have become aware that the recent paper 
of Bordag, Kirsten, and Vassilevich arrives at results closely related 

The authors of 



to those obtained by Barton [15 



241 make use of field 



theoretical path integral methods which are in themselves quite different 



from those used in |15[ as well as in the present paper, but the expression 



for the Casimir energy obtained in Section IV in |24| nevertheless parallels 
that obtained in ||15[| . (There are some differences in numerical coefficients of 
the divergent terms, due to different regularization methods employed.) We 
thank the authors of [2~1| for making us aware of this correspondence. 
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